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BEST UNIFORM RATIONAL APPROXIMATION OF x a ON [0,1] 



Herbert Stahl 

Abstract. A strong error estimate for the uniform rational approximation of x a 
on [0, 1] is given, and its proof is sketched. Let E nn (x a , [0, 1]) denote the minimal 
approximation error in the uniform norm. Then it is shown that 



lim e 2n ^ an E nn (x a ,{0,l]) = 4 1+Q | shi7ra| 

n — >oo 

holds true for each a > 0. 

1. Introduction 

Let II n denote the set of all polynomials of degree at most n £ N with real 
coefficients; TZ mn denote the set {p/q\p £ II m , q £ II n , q ^ 0}, m,n £ N, of 
rational functions; and the best rational approximant r* m £ TZ m n, m,n £ N, and 
the minimal approximation error E mn — E mn (x a , [0, 1]) be defined by 

(1-1) E mn := \\x a - r* mn \\ [0A] = inf \\x a - r|| [0 ,i], 

where || • ||k denotes the sup norm on K C M. It is well known that the best 
approximant r* mn exists and is unique within lZ mn (cf. [Me, §§9.1, 9.2] or [Ri, 
§5.1]). The unique existence also holds in the special case (n = 0) of best polynomial 
approximants. 

Since f a {x) := \x\ a is an even function on [—1, 1], the same is true for its unique 
approximant r* m = r* nn (f a , [— 1, 1]; •), and consequently a substitution of z 2 by 
z shows that approximating \x\ 2a on [—1,1] and x a on [0,1] poses an identical 
problem. We have 

(1.2) £ 2m ,2„(M 2Q ,[-l,l]) = E mn (x a : [0,l}) forallm,neN. 

From Jackson's and Bernstein's theorems about the interdependence of approx- 
imation speed and smoothness of the function to be approximated (cf. [Me, §§5.5, 
5.6]) we know that in case of a € K.+\N the minimal error E m fi{\x\ a , [—1, 1]) be- 
haves like 0(m~ a ) as m — > oo. In [Bel, Be2] Bernstein proved a result that is 
deeper and much more difficult to obtain; he showed that the limit 

(1.3) lim m a E mfi (\x\ a ,[-l,l]):=f3(a) 
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exists for each a > 0; however, an explicit expression for the constant (3(a), a > 0, 
still is not known. 

In case a = 1 the number (3 := [3(1) — 0.28016... is known as Bernstein's 
constant. In [Bel, p. 56] Bernstein raised the question whether (3 can be expressed 
by known transcendentals or whether it defines a new one, and based on numerical 
upper and lower bounds for (3, which he calculated up to a precision of ±0.005, he 

made the tentative conjecture (3 = 1/(2-^), which, however, has been disproved in 
[VC1] by extensive and nontrivial high precision calculations. 

For large values of a, Bernstein was able to establish in [Be2] an asymptotic 
expression. He showed that 

(3(a) 1 
hm — — = — . 

a^oo T(a)\ sin(7ra/2)| n 

While Bernstein's investigations on best polynomial approximation of |a;| and 
\x\ a were published in the period between 1909 and 1938, the study of best rational 
approximation of \x\ a was only started in 1964 by Newman's surprising (at the 
time) result in [Ne] that 

(1.5) ±e- 9 ^ < E nn (\x\,[-l,l]) <3e-^ for all n = 4, 5,.... 

A comparison of this result with (1.3) shows that the convergence behavior of 
rational approximants is essentially better than that of polynomials. Ncwmann's 
investigation has triggered a whole series of contributions, from which we select a 
short list with papers that contain substantial improvements of the error estimate 
in the uniform norm. 



E nn (x a , [0, 1]) < e - c (")^\ a e R+ in [FrSz]; 

E nn (x^ 3 ,[0,l}) <e- c ^, in [Bui]; 

E nn (x a , [0, 1]) <e- c(a ^ 7 a G R+, in [Gol]; 

l e -^ < En^X 1 / 2 , [0, 1]) < e -V^(l-0(n-^)) t j n [Bu2] . 

e -c(a)v^ < Enn ( x « j [ 0j 1]), a e Q + \ N) i n [Go2]; 

^v^a+e) < E nn (x a , [o,i]) < e - 7r ^( 1 - £ ), 

a e K+\N, e > 0, n > n {e, a) in [Go3]; 

E nn (x^ 2 , [0,1]) <cne-^, in [Vjl]; 

i e -*V*i < E nn (x x '\ [0, 1]) <ce-*^ in [Vj2]; 

e -d(«)VH < E nn (tfx, [0, 1]) < e - c ^ s ^, s e N in [Tz]. 



Here c, c(a), . . . denote constants. Relation (1.2) allows us to transfer these results 
to the problem of approximating on [—1,1]. 

It may be appropriate to repeat a remark from [Go2], where it was pointed out 
that Newman's result can be obtained rather immediately from an old result (from 
1877) by Zolotarev. In this sense, the investigation of rational approximants dates 
back even further than Bernstein's work on polynomial approximants. 



(1.4) 
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The best result known so far for the rational approximation of x a on [0, 1] was 
obtained independently by Ganelius [Ga] in 1979 and by Vjacheslavov [Vj3] in 1980. 
They proved that for a £ R + \N there exists a constant c\ = c\(a) > such that 

(1.6) liminfe 2 ^\E ml (x tt ,[0,l]) > ci(a) 

n^oo 

and, conversely, that for each positive rational number a e Q + there exists a 
constant c 2 = 02(a) < 00 such that 

(1.7) \imsupe 27r ^E nn (x a ,[0,l]) < C2 (a). 

n^oo 

Both authors were not able to show that C2 = 02(0;) depends continuously on a. 
Thus, inequality (1.7) remained open for a e K+\Q; however, in [Ga] Ganelius was 
able to prove the somewhat weaker estimate 

(1.8) E nn (x a ,[0,l})<c 2 (a)e 2 ^ +c ^^ for n > n (o 2 (a), 03(a)), 

which holds for each a > 0. In (1.8) 02(a) and 03(a) are constants depending on a. 

The results (1.6)— (1.8) give the correct exponent —2ny/an in the error formula; 
however, nearly nothing is said about the coefficient in front of the exponential term. 
The determination of this coefficient is the subject of the present note. Practically as 
a byproduct, we prove the upper estimate (1.7) for irrational exponents a € R+\Q. 

2. The result 

Theorem 1. The limit 

(2.1) lim e 27r ^E nn (x a , [0, 1]) = 4 1+Q | sin7ra| 

n — >oo 

holds for each a > 0. 

Remarks. (1) From (2.1) we deduce that the approximation error E nn (x a , [0,1]) has 
the asymptotic behavior 

(2.2) E nn (x a ,[0,l}) = 4 1+a | sin7ra|e- 27r ^"(l + o(l)) asn^oo, 
and equivalently it follows with (1.2) that 

(2.3) E nn (\x\ a ,[-l,l\) =4 1+Q / 2 |sin7ra/2|e-' r ^ : (l + o(l)) asn^oo 
for each a > 0. 

(2) Not only the explicit expression on the right-hand side of (2.1) but already the 
existence of the limit represents a result difficult to obtain. The value 4 1+Q | sin7ra|, 
a > 0, is the analogue of Bernstein's constant (3(a) in (1.3) for the case of ratio- 
nal approximation. It has already been noted in the introduction that an explicit 
expression for /3(a) is still not known. The best we know is Bernstein's asymp- 
totic formula (1.4). Since in (1.4) we have considered approximation on [—1,1], 
the counterpart of the asymptotic value lr(a)| sin7ra/2| for /3(a) is the value 
41+0/2 1 3^(^0-/2)! in case of rational approximation. 
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(3) If we turn our attention to the special case of rational approximation of \x\ 
on [—1, 1], then it follows from (1.2) that 

(2.4) £ 2n , 2 „(| a; |,[-l,l])=-B„n(\/i,[0,l]) forneN, 
and hence we deduce from (2.1) that 

(2.5) lim e-*^E nn (\x\,[-l,l]) = 4 1+1 / 2 | sin f | = 8. 

n— >oo 

Limit (2.5) has been conjectured in [VRC] on the basis of high precision calculations 
and was proved in [St] . 

It may be surprising that in case of rational approximation, which is in many 
respects more complex than the polynomial case, limit (2.5) has a rational value, 
while in the polynomial case Bernstein's question in [Bel] about the character of 
the number (3 — (3(1) is still open and the numerical results in [VC1] show that (3 
cannot be a rational number with a moderately small denominator. 

3. Outline of the proof of Theorem 1 
If limit (2.1) is proved for one of the paradiagonal sequences 

{JW(* a >[o,i])}~| fc |, 

fceZ fixed, then it holds also for the diagonal sequence 

[o,i])}~ =1 . 

It turns out that 

(3.1) n + k = m:=n+l + [a], aeR+\N,neN, 

is a good choice for the numerator degree m. From the theory of best rational 
approximants we learn that the error function 

(3.2) e n (z):=z a -r* mn (z), z e C\R_, 

has exactly 2n + 2 + [a] zeros in the interval (0, 1). Hence the theory of multipoint 
Pade approximants is applicable, and it gives us rather precise information about 
the structure of the numerators and denominators of r* mn (cf. [GoLa; StTo, §§6.1, 
6.2]). 

In the next step the error function e„ and the approximant r* mn will be trans- 
formed in such a way that the resulting function \I>„ is analytic in C\R with possible 
exceptions in a disc A(i?) with radius R > around the origin. The function \P n 
has boundary values from both sides of M\A(i?) that allow a comparison with a 
special logarithmic potential. The potential will be introduced after Theorem 2 
below. 
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The transformation of e„ and r* mn into the function \I>„ is carried out in several 
steps. The intermediate functions are defined as 

(3.3) rJz) := ^£Z ; 

Aw 2a — 1 1 

(3.4) R n (w) := ——r n (el/ a w) - — , e n := E mn (x a , [0, 1]); 

(3.5) $„H := ^ (i?„H + v^nH 2 - 4) ; 

i/}($> n (w)) for Im(w;) > 0, 



(3.6) *„(«;) := 



with 



V>($„(w;)) for Im(w;) < 0, 



(3.7) V(^) := 7 r - 

sin not — i(cosira)z 

In (3.4) a new variable w is introduced implicitly by 

(3.8) w:=s- 1/a z, zeC. 

The properties of each new function r„, R n , $„ have to be studied carefully. The 
properties of the last function ^>„ is summcrized in 



Lemma 1. The function ^> n is analytic and different from zero in C\(KU A(R)), 
R > appropriately chosen, and there exist constants c\ , . . . , C4 such that 

(3.9) |log|*„H|| < Cl \w\~ 2a forweM.-\A(R),n>no( Cl ,R); 

(3.10) |log|*„(w)w a 4sin7ra|| < c 2 \w\~ a for w e R+\A(R),n > n Q (c 2 ,R); 

(3.11) |log|*„H|| < c 3 forwedA(R),n>n a {c 3 ,R). 

If we consider the representation 

(3.12) log |* n («;) I =rl> n (w)- J g D (w,t)dii n (t), 

where D := C\(R_ U A(R)), go(w,t) the Green function in D, ip n a harmonic 
function in D with 

(3.13) VnH = log|*(to)| forwedD, 
and \i n a measure with 

(3.14) supp(/i„) C [R, 00] and ^ n >0 on [R, £~ 1//o! ], 
then in addition to (3.9)-(3.11) we have the inequalities 

(3-15) IKI [e -i/«.,oo]ll<l. 

(3.16) \fi n ([R,e~ 1/a ])-2n\<c 4 for n > n (c 4 , R). 

Remark. Estimates (3.9), (3.10), and (3.16) contain the information that is most 
relevant for the proof of Theorem 1 . 

The function log |*„| will be compared with a special Green potential p n . The 
definition of this potential is based on the following. 
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Theorem 2. For each a E (l,oo) there uniquely exists a Green potential 
(3.17) Pa(w):= j g D (w,t)du a (t), D :-- 



with 
(3.18) 

that satisfies 
(3.19) 



v a > 0, supp(^ a ) = [b a , a], < b a < a, 



Pa(w)| 



= logw for w e [b a ,a], 
> \ogw for w € (0,6 a ). 



For the constant b a and the measure v a appearing in (3.17)-(3.19) we have 

(3.20) b a -» V2 
and 

(3.21) asa^oo. 

A proof of Theorem 2 can be derived from [St, Theorem 2]. It is necessary to 
change the domain of definition D by the transformation w i— > 1/ -^/w. Basic tools 
in the proof of Theorem 2 are estimates for certain elliptical integrals. 

The potential p n for a comparison with the function log is now defined as 



(3.22) 
with 

(3.23) 



p n (w) := -ap a (cw) 



a := 



■ sin ira 



l/a 



and c := |4sin7ra 



1/q 



From (3.21), (3.23), and (3.16) together with other information provided by Lemma 
1 it then follows that 



(3.24) 



4 sin7ra 



l/a 



;7 r 1 /2(2n)/ Ct 



4 sin ira\ 



as n — > oo. 



Theorem 1 then immediately follows from (3.24). 

What we have given here is only a sketch of the overall structure of the proof of 
Theorem 1. Some of the steps demand rather subtle analysis. 

Note Added in Proof. In [VC2] Varga and Carpenter have calculated numerical 
approximations for the right-hand side of (2.1) for the six values a = j/8, j = 
1,2,3,5,6,7, and conjectured formula (2.1) on the basis of these numerical val- 
ues. The numerical results for the two cases a = 1/4 and a = 3/4 are especially 
interesting since here the right-hand side of (2.1) is rational. 
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